Shape modeling technique KOALA validated by ESA Rosetta at (21) Lutetia 
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J — Abstract 

We present here a comparison of our results from ground-based observations of asteroid (21) Lutetia with imaging data acquired 
I— I during the flyby of the asteroid by the ESA Rosetta mission. This flyby provided a unique opportunity to evaluate and calibrate 
our method of determination of size, 3-D shape, and spin of an asteroid from ground-based observations. Knowledge of certain 
observable physical properties of small bodies {e.g., size, spin, 3-D shape, and density) have far-reaching implications in furthering 
(— I our understanding of these objects, such as composition, internal structure, and the effects of non-gravitational forces. We review 
CL(the different observing techniques used to determine the above physical properties of asteroids and present our 3-D shape-modeling 
Q technique KOALA - Knitted Occultation, Adaptive-optics, and Lightcurve Analysis - which is based on multi-dataset inversion. 
^ We compare the results we obtained with KOALA, prior to the flyby, on asteroid (21) Lutetia with the high-spatial resolution 
images of the asteroid taken with the OSIRIS camera on-board the ESA Rosetta spacecraft, during its encounter with Lutetia on 
2010 July 10. The spin axis determined with KOALA was found to be accurate to within two degrees, while the KOALA diameter 
determinations were within 2% of the Rosetta-derived values. The 3-D shape of the KOALA model is also confirmed by the 
spectacular visual agreement between both 3-D shape models (KOALA pre- and OSIRIS post-flyby). We found a typical deviation 
of only 2 km at local scales between the profiles from KOALA predictions and OSIRIS images, resulting in a volume uncertainty 
provided by KOALA better than 10%. Radiometric techniques for the interpretation of thermal infrared data also benefit greatly 
from the KOALA shape model: the absolute size and geometric albedo can be derived with high accuracy, and thermal properties, 
for example the thermal inertia, can be determined unambiguously. The corresponding Lutetia analysis leads to a geometric albedo 
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of 0.19 + 0.01 and a thermal inertia below 40 J m s"" K" , both in excellent agreement with the Rosetta findings. We consider 
this to be a validation of the KOALA method. Because space exploration will remain limited to only a few objects, KOALA stands 
as a powerful technique to study a much larger set of small bodies using Earth-based observations. 
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■1. Remote-sensing shape modeling 

Perhaps the most striking observable of any asteroid is its 
shape. In 1993, spacecraft exploration revealed for the first 
time the stunning non-spherical shape of asteroid (951) Gaspra 
when NASA's Galileo spacecraft made the first of its two aster- 
oid encounters, on its way to Jupiter. Asteroids had remained 
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point-sources in the sky since the discovery of (1) Ceres in 
1801 by Piazzi, almost two centuries before. Only the ad- 
vent of space exploration and large Earth-based telescopes (e.g., 
Arecibo and Goldstone radio telescopes, space-based-optical 
- HST, or ground-based near-IR, equipped with adaptive op- 
tics - Palomar, Lick, CFHT, Keck, ESO VLT, Gemini) allowed 
their apparent disks to be spatially resolved, and their irregu- 
lar shapes to be imaged. The past decade has seen a revolution 
in our understanding of the physical properties {e.g., size, 3-D 
shape, spin axis) of asteroids. This revolution has come about 
thanks to improved observing facilities and, of equal impor- 
tance, from improved methods of analysis. 

Determination of the physical properties for a statistically 
relevant set of asteroids has many implications for our under- 
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standing of these remnants of solar-system formation and, in 
turn, can be expected to improve our understanding of the his- 
tory and evolution of the Solar System. For instance, the distri- 
bution of spin axes of the larger asteroids (diameter larger than 
~100 km) on the celestial sphere is not expected to be isotropic. 
Numerical hydrocode simulations have predicted a slight ex- 
cess in prograde rotators, due to the gas-pebble interaction in 
the protoplanetary disk (Johansen and Lacerda 2010). Simi- 
larly, the spin state of small asteroids (diameter not larger than 
few kilometers) is dominated by the non-gravitational YORP 
effect (Hanus et al. 201 1). Statistical knowledge of spin coordi- 
nates, how they are distributed within and among asteroid fam- 
ilies of different ages, will provide strong constraints on the ef- 
fectiveness of YORP (Slivan et al. 2003). 

Reconstruction of the 3-D shape (including the size) is re- 
quired to estimate the volume of an asteroid, which in turn is 
used to derive its density, possibly the property most funda- 
mental to our understanding of an asteroid (Britt et al. 2002). 
Observations of the surface of an asteroid, such as colors, spec- 
tra, or phase effects, can provide clues to the surface compo- 
sition of the asteroid. This may or may not be related to the 
bulk composition of the body (e.g., Elkins-Tanton et al. 201 1). 
Masses for asteroids can be determined from a spacecraft flyby, 
from the orbital motion of a natural moon, or even from the 
perturbations of asteroids on other bodies, such as Mars (Hilton 
2002). In most cases, however, the uncertainty in the density 
is dominated by the uncertainty in the volume, rather than the 
uncertainty in the mass (Merline et al. 2002). Precise recon- 
struction of the 3-D shape is therefore of high importance for 
all asteroids for which a mass has been, or will be, estimated 
(e.g., Hilton 2002; Mouret et al. 2007; Baer et al. 2011). 

From the comparison of an asteroid's density with the densi- 
ties of its most-likely constituents, we can constrain the macro- 
porosity (large-scale voids) in its interior, probably produced 
by impacts over its history (Britt et al. 2002). These impacts 
could have partially disrupted the body, producing large-scale 
fractures, or even totally disrupted the body, with subsequent 
reaccumulation of the resulting fragments, leading to a "rubble- 
pile" structure. 

Evidence of gigantic, but less than totally disruptive, im- 
pacts can be seen by high-resolution imaging and also in- 
ferred from our shape-modeling of asteroids. The huge impact 
craters evident in the images of C-type asteroid (253) Mathilde 
(Veverka et al. 1997) are thought to be about as large as could 
be sustainable by a body without disruption. There have been 
suggestions that these craters were created by compaction of 
low-density target material, rather than explosive ejection typ- 
ical of hard-rock impacts (Housen et al. 1999). Already, our 
ground-based adaptive-optics imaging has shown what appear 
to be facets or depression-like features, similar to those seen 
on Mathilde, in some other large C-type asteroids (e.g., (511) 
Davida in Conrad et al. 2007). Alternatively, some of our other 
images of C-type asteroids, such as (52) Europa, appear to bear 
no evidence of giant impacts (see Merline et al. 2011). Evalu- 
ation of the prevalence of such large impact events can give us 
insight into the size and frequency of these impact events over 
time, and thus into the history of the impacting population. 



We summarize below the most common of the many observ- 
ing techniques used to derive size, 3-D shape, and spin-vector 
coordinates and highlight some of their advantages and draw- 
backs. Then, in Section 2, we describe our KOALA multi-data 
shape-modeling algorithm. In Section 3, we present a com- 
parison of the results produced by KOALA (from Earth-based 
observations) with those derived from the ESA Rosetta flyby of 
asteroid (21) Lutetia. In Section 4, we use our KOALA model 
in conjunction with mid-infrared data and a thermophysical 
model to derive the thermal properties of Lutetia and compare 
the results with those derived using thermal observations from 
the Rosetta spacecraft (Gulkis et al. 201 1) and from the ground, 
making use of the shape model from the flyby (O'Rourke et al. 
201 1). We assess the accuracy of the KOALA shape-modeling 
method in Section 5. 

1.1. Optical lightcurve 

Historically, spin properties and triaxial-ellipsoid shapes 
have been studied largely through observations of rotationally 
induced variability in disk-integrated brightness (lightcurves). 
Indeed, the object's shape, its rotational state (period and spin- 
vector coordinates), and the scattering properties of its surface 
can be determined from the analysis of its lightcurves over time 
(as the viewing/illumination geometry changes). 

For about a decade, starting with the lightcurve inversion 
algorithm presented by Kaasalainen and Torppa (2001) and 
Kaasalainen et al. (2001), lightcurves also have been used ex- 
tensively to derive 3-D shape models of asteroids (see the ex- 
amples for asteroids Gaspra and Steins in Fig. 1). The 3-D 
shape models and spin properties of more than 200 asteroids 
already have been derived (these are accessible from DAMIT ^, 
see Durech et al. 2010). These shape models are, however, lim- 
ited to dimensionless, convex shapes, with limited spatial reso- 
lution. Recently, Durech et al. (201 1) have shown that the size 
of these models can be set by using other types of input data 
(e.g., stellar occultation profiles, see Sect. 1.2). The intrinsic 
convex nature of the models precludes accurate determination 
of the volume, and hence density, of the objects, however 

Because lightcurve observations require neither large tele- 
scope aperture nor specialized instrumentation, they are, and 
will remain, a major source of information on small bodies. 
Thousands of lightcurves, for hundreds of asteroids^, have 
been accumulated during the last half century. Amateur as- 
tronomers contribute significantly to this ever-growing database 
of lightcurves {e.g., Behrend et al. 2006; Durech et al. 2007). 

Sparse photometry (i.e., when the typical separation between 
measurements is larger than the rotation period, as opposed to 
historical lightcurves, which are dense in time) can also be used 



-DAMIT: http://astro.troja.mif.cuni.cz/projects/asteroids3D/web.php 
^e.g., the Asteroid Photometry Catalogue (APC) or the Asteroid Lightcurve 

Database (LCDB, Warner et al. 2009) list more than 6000 lightcurves for about 

700 asteroids: 

APC: http://asteroid.astro.helsinki.fi/apc 

LCDB: http://www.minorplanet.info/Hghtcurvedatabase.html 

about 2300 lightcurves for more than 1700 asteroids have been ac- 
quired by the CdR group: http://obswww.unige.ch/~behrend/page_cou.html 
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Figure 1: Comparison of four sliape models, derived using ground-based ob- 
servations, with high-spatial-resolution images obtained in situ by spacecraft 
(left-to-right): (951) Gaspra lightcurve inversion model from Kaasalainen et al. 
(2()()2a), image from NASA Galileo; (2867) Steins lightcurve inversion model 
from Lamy et al. (2008a, b), image from ESA Rosetta; (25143) Itokawa radar 
delay-Doppler model from Ostro et al. (2005), image from JAXA Hayabusa; 
and (21) Lutetia KOALA model from CaiTy et al. (2010b) and Drummond et al. 

(2010) , image from ESA Rosetta. 

to reconstruct 3-D models (see Kaasalainen 2004). Hanus et al. 

(201 1) have used a combination of sparse photometry, together 
with dense lightcurves, to derive about 100 new shape models, 
using measurements extracted from large all-sky surveys (such 
as USNO, Catalina, Siding Spring, and Hipparcos). Knowledge 
of the absolute photometry (as opposed to relative photometric 
measurements, as in a dense lightcurve) is, however, required to 
make these sparse measurements useful to 3-D shape modeling 
(see Hanus et al. 2011, for a detailed discussion). From the 
upcoming PanSTARRS and Gaia surveys we can expect hun- 
dreds of thousands of objects to be modeled using this method 
(Durech et al. 2005). 

1.2. Stellar occultation 

Occultations occur when a solar-system object passes be- 
tween a star and Earth, causing the asteroid's shadow to cross 
some portion of Earth as a track. Because the star is very far 
relative to the Earth-asteroid separation, the shadow cast by the 
asteroid is effectively parallel. Thus, the width of the shadow 
track (perpendicular to the track) gives the maximum width of 
the asteroid in the cross-track direction. It is usually not possi- 
ble to get a high density of observers stationed across the track, 
and thus this dimension may not be so well established. But 
along the track, the size of the asteroid, as an along-track chord 
on the asteroid, corresponding to the position of the observer 
within the shadow, is given by the length of time of the bhnk- 
out event. With many observers, many along-track events can 
be recorded and the blink-out intervals are converted to chord 
lengths at the asteroid by using the asteroid's known speed (see 
Millis and Dunham 1989, for a review). 

One advantage of stellar occultations is that very small mi- 
nor planets can be probed (provided one accounts properly for 
diffraction effects). Even stellar occultations by small trans- 
Neptunian objects (TNOs) of a few kilometers diameter can be 
recorded (see Roques et al. 2009). Finally, stellar occultations 
provide a powerful means for studying thin atmospheres and/or 
exospheres {e.g., Sicardy et al. 2003). 

In theory, three occultation events (each recorded by enough 
observers to secure a unique apparent-ellipse to be fit on the 



chords) provide enough constraints to determine the triaxial di- 
mensions (ellipsoid) and spin-vector coordinates of an asteroid 
(Drummond and Cocke 1989; Dunham et al. 1990). The num- 
ber of chords that result from an event will often be larger with 
brighter occulted stars. This is because many more observers 
can be fielded if the required telescope aperture is modest. 
In practice, however, occultations of bright stars by any given 
asteroid occur rather infrequently. This difficulty in obtaining 
sufficient chords, plus the noise level often encountered (espe- 
cially the systematic errors coming from imperfect knowledge 
of asteroid and star astrometry, combined with observer timing 
errors), strongly limit our ability to construct 3-D shape mod- 
els and derive spin properties from occultations alone. Never- 
theless, stellar occultations are an efficient way to provide ad- 
ditional size/shape information, particularly for shape models 
that otherwise lack a scale, such as those from lightcurves alone 
(see 1.1 and Durech et al. 201 1). 

With the availability of low-cost GPS positioning equipment 
and CCD cameras (the majority of measurements are made by 
amateur astronomers), the accuracy of occultation timings has 
improved greatly over the last decade. From this improved pre- 
cision, together with the refinement of the orbits of small bod- 
ies expected to result from Gaia/PanSTARRS (allowing an ex- 
tremely precise prediction of the occultation track on Earth, see 
Tanga and Delbo 2007), we can expect that stellar occultations 
will become ever more useful in the determination of certain 
physical properties of asteroids, especially for objects having 
small angular diameters. 

1.3. Thermal radiometry 

The amount of thermal emission from an asteroid is mainly 
a function of its physical diameter and surface albedo, and, to a 
lesser extent, the physical properties of its surface {e.g., thermal 
inertia, roughness). Main-belt asteroids are among the bright- 
est sources in the sky in the mid-infrared (5-20 microns), so in- 
frared satellites (IRAS, ISO, AKARI, Spitzer, Herschel, WISE) 
have been able to acquire observations of a vast number of these 
objects (see Mainzer et al. 201 1, for instance). 

Estimates for size and albedo are model-dependent, how- 
ever, and several thermal models co-exist, from the simple 
Standard Thermal Model (STM) of non-rotating spheres of 
Lebofsky et al. (1986) to the detailed Thermophysical Model 
(TPM) of Lagerros (1996, 1997), having a complete descrip- 
tion of 3-D shape and surface properties. Unfortunately, the 
systematic uncertainties involved when applying those models 
(resulting from their respective assumptions and approxima- 
tions) are not always properly taken into account in estimating 
error bars, and results often differ from one determination to 
another by more than the quoted uncertainties (see Table 3 in 
Delbo and Tanga 2009, illustrating the issue). For instance, it 
now seems that the database of 2 228 diameters (Tedesco et al. 
2002) estimated using the overly simple STM with IRAS data 
may be biased by a few percents (see the re-analysis of IRAS 
data by Ryan and Woodward 2010). 

The radiometric technique to derive sizes and albedos from 
thermal infrared data benefits greatly from the availability 
of 3-D shape models. For example, the absolute size of 
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(25143) Itokawa derived by Mulleretal. (2005), based on a 
3-D shape model implementation in a TPM code, agreed to 
within 2% of the final in-situ result from the Hayabusa mis- 
sion (Fujiwara et al. 2006). In general, mid-infrared observa- 
tions are highly valuable for scaling dimensionless shape mod- 
els (similar to the situation with stellar occultations, as men- 
tioned in Sect. 1 .2) and mid-infrared data are available for sev- 
eral thousand asteroids. These data even allow determination 
of the most likely spin-axis solutions in cases where lightcurve 
inversion techniques lead to more than one possible shape and 
spin-axis solution (see, e.g., Miiller et al. 2011, for a recent ex- 
ample). In cases where the shape-model already comes with 
size information (or alternatively, if many thermal observations 
are available for a given target), it is possible to derive the ther- 
mal inertia, indicative of the surface characteristics: e.g., bare 
rock, ice, boulders, dust regolith (see, for instance Miiller et al. 
2005; Delbo et al. 2007, among many others). 

In some particular cases of extensively observed asteroids, 
the thermal radiometry can also provide hints on the 3-D shape, 
through the measure of the apparent projected cross-section of 
the asteroid on the plane of the sky at each epoch. By compar- 
ing the predicted with observed thermal fluxes of Lutetia un- 
der many geometries, O'Rourke et al. (2011) have shown that 
adding a hill/plateau, whose size remains within the quoted 3-D 
shape uncertainty, could explain the discrepancies observed for 
a certain observing geometry. 

1.4. Radar delay-Doppler echoes 

Radar observations consist in transmitting a radio signal to- 
ward the target and recording the echo. The distribution of the 
echo power in delay time and Doppler frequency is related to 
the spin state and 3-D shape of the target (see the reviews by 
Ostro 1989; Ostro et al. 2002). The time and frequency preci- 
sion currently achievable (Arecibo, Goldstone) allow the study 
of very small objects, the main limit of radar observations be- 
ing the distance of the target (echo power scales inversely with 
distance to the fourth power). This is why most radar studies 
have concentrated on Near-Earth Objects (NEOs), while dedi- 
cated observations of Main-Belt Asteroids (MBAs) have been 
more hmited (see Ostro et al. 2002). 

The difficulty in constructing 3-D shapes from a series of 
delay-Doppler "images" is due to the absence of a direct, one- 
to-one, link between each region of the surface with a pixel in 
delay-Doppler space. Indeed, all points situated at the same 
range from the observer, and moving at the same radial velocity 
(possibly spread over the apparent disk) will contribute to a sin- 
gle delay-Doppler pixel. So delay-Doppler images are many- 
to-one maps of the shape as highlighted by Ostro et al. (2002): 
there is no a priori way to determine how many regions will 
contribute to a single pixel, which strongly contrasts with the 
one-to-one mapping ("WYSIWYG") achieved in disk-resolved 
imaging. 

Radar echoes remain, however, the best way to determine the 
physical properties of NEOs {e.g., the very small NEO Itokawa 
in Fig. 1). For instance, the possible migration of the regolith at 
the surface of a fast-rotating asteroid triggered by YORP spin- 



up (Walsh et al. 2008) was suggested by the detailed 3-D shape 
of the NEO (66391) 1999 KW4 (Ostro et al. 2006). 

7.5. Disk-resolved imaging 

Since the 1990s, with the advent of the Hubble Space 
Telescope and large ground-based telescopes equipped with 
adaptive optics (AO: Lick, CFHT, Keck, ESO VLT, Gemini), 
we have access to the angular resolution required to resolve 
the apparent disk of asteroids {e.g., Saint-Pe et al. 1993a,b; 
Drummond et al. 1998; Parker et al. 2002, 2006). 

From a time-series of disk-resolved images, spin-vector 
coordinates can be derived (using previous knowledge of the 
rotation period) by analyzing the changes in the apparent shape 
of the asteroid profile (see for instance Thomas et al. 1997b; 
Drummond and Christou 2008), or by following the apparent 
path taken by an albedo patch on the surface during the rotation 
{e.g., Thomas et al. 2005; Carry et al. 2008). Triaxial shapes 
(ellipsoids) can also be derived (see Drummond et al. 2009; 
Schmidt et al. 2009, for instance), and topography (such as the 
presence of facets or craters) studied from profile deviations 
to the ellipsoid {e.g., Thomas et al. 1997a; Conrad et al. 2007). 
With sufficient spatial resolution, imaging of the asteroid disk 
can allow construction of albedo maps of the surface, allowing 
the study of composition heterogeneity or localized space 
weathering effects {e.g., Binzel et al. 1997; Li et al. 2006, 
2010; Carry et al. 2008, 2010a). In the case of asteroids visited 
by spacecraft, high spatial-resolution and precise photometry 
can be used to derive precise shape and digital terrain models 
by using stereophotoclinometric techniques (see examples in 
Gaskell et al. 2008). 

The size and 3-D shape resulting from disk-resolved images 
are of great value, being obtained in a direct manner (as com- 
pared to an indirect determination of the shape from lightcurve 
inversion, for instance). The best angular resolution ^ of current 
Earth-based telescopes is about 0.04". Due to systematics, 
however, we have found that our ability to accurately measure 
sizes degrades rapidly below about 0. 10", based on simulations 
and observations of the moons of Saturn and other objects 
(Carry 2009; Drummond et al. 2009). The sample of asteroids 
observable {i.e., having angular sizes that get above about 
0.10") is therefore limited to about 200 (over a given 10-year 
span). 



1.6. Interferometry 

Apart from building larger telescopes, one efficient way to 
improve the angular resolution is to combine light beams from 
separated telescopes and to observe the resulting interference 
(fringes). In such a mode, each telescope aperture plays the role 
of a sample aperture within a virtual telescope whose extent is 
the largest distance between the two telescopes (the spatial res- 
olution is still given by the equation in footnote 5, except that 



^limited by the diffraction, wliicli acts as a low-pass filter with a cutoff fre- 
quency approximated by = AjD (radian), A being the wavelength and T) the 
diameter of the telescope aperture. 
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D is now the distance between the apertures). For instance, 
with telescopes separated by about 80 m, the VLTI provides an 
order-of-magnitude improvement in angular resolution over a 
single telescope of the VLT (8 m aperture). 

This improvement in the resolution, however, occurs at the 
price of a loss in complete spatial information, because the 
virtual aperture is under-sampled. The very high angular- 
resolution is limited, at a given instant, to a single line on 
the plane of the sky, which is given by the baseline linking the 
two apertures. Along that baseline, the signal is directly related 
to the Fourier transform of the flux distribution on the plane 
of the sky. To expand the coverage in the spatial-frequency 
domain, and thus to allow the construction of 2-D images of 
the target, one must increase the number of projected baselines. 
This is commonly achieved by either adding multiple physical 
baselines (i.e., adding more telescopes) and/or by making ob- 
servations throughout the night, as the Earth's rotation causes 
a progression of the position angle of the baseline on the sky. 
However, because asteroids typically complete their rotation in 
only a few hours, the method using multiple telescopes is more 
effective at providing higher spatial sampling. Below we de- 
scribe another promising technique to effectively increase the 
number of baselines, by using a Fizeau design. 

For objects with no a priori information, the overall dimen- 
sions and spin properties can be determined under the assump- 
tion that the shape is well-described by a triaxial ellispoid 
(Li et al. 2011). For example, the two orthogonal directions 
of the Fine Guidance Sensor (FGS) on-board HST have been 
used to study binarity and measure the size of several MBAs 
(Tangaetal. 2001; Hestroff'er et al. 2002; Tangaetal. 2003). 
The fringes of interference, however, also contain informa- 
tion on the apparent shape, and can be used in combination 
with other data to derive 3-D shape models. For instance, 
Kaasalainen and Lamberg (2006) used FGS data in combina- 
tion with optical lightcurves to refine the 3-D shape model of 
(15) Eunomia from Kaasalainen et al. (2002b). 

Interferometry in the mid-infrared has been used recently 
also, with promising results. In this wavelength regime, the 
signal is linked to the distribution of temperature (i.e., emit- 
ted light) on the surface of the target (as opposed to the re- 
flected light seen at visible wavelengths). Delbo et al. (2006) 
and Delbo et al. (2009) have combined thermal infrared TPM 
models (see Sect. 1.3) with interferometry, allowing the size 
of 3-D shape models to be set (similar to stellar occupation and 
thermal-infrared-only, see Sections 1.2 and 1.3 and Matter et al. 
2011). 

The current limitation of interferometry is set by the sensitiv- 
ity of available facilities and is driven largely by the integration 
time, which is usually short. Because the light beams from the 
two apertures traverse different paths in the atmosphere, their 
wavefronts encounter different turbulence-dominated perturba- 
tions, leading to a shift in their phase. Delay lines are used to 
"slow down" one beam with respect to the other, and to match 
their phase. The technical difficulties in maintaining coherence 
in this process limits the integration time to few hundredths or 
thousandths of a second. A new generation of instruments (like 
PRIMA at VLT) with fringe-tracking systems will overcome 



this limitation in the near future, allowing fainter sources to be 
targeted. 

An interesting compromise between co-axial interferome- 
ters and traditional, filled aperture, telescopes can be found in 
Fizeau-type instruments such as LINC-NIRVANA, being built 
for the Large Binocular Telescope (LBT, see Hill et al. 2010), in 
which improved coverage of the spatial-frequency plane (com- 
monly called uv-plane, with u and v standing for the two or- 
thogonal unit directions) can be achieved in less telescope time. 
This design (Fizeau vs. pupil-plane interferometer like the 
VLTI) allows instantaneous filUng of the uv-plane up to the fre- 
quency set by the 8 m apertures (Fizeau interferometers are true 
imaging devices and produce direct images of the plane of the 
sky), with an additional filling along one dimension up to the 
frequencies corresponding to the maximal baseline (22.7 m for 
LBT). Two or three epochs, with different position angles of the 
baseline on the plane of the sky, will be enough to fill the uv- 
plane up to the maximal extent of the telescope (achievable for 
transiting asteroids, especially with high elevation, when the 
position angle evolves quickly). For this purpose, therefore, 
the LBT is equivalent to a 22.7 m telescope with a mask, cor- 
responding to the configuration of component apertures, placed 
in the entrance pupil. 

The opening of interferometric studies to longer wavelengths 
also has great potential, with the use of millimeter and sub- 
millimeter arrays, where there are fewer technological limita- 
tions than in the optical range. Future facilities such as ALMA, 
with 50 antennas (translating into 1 225 baselines vs. only 
2baseHnes for MIDI at the VLTI, see Delbo et al. 2009), wiU al- 
low dense spatial coverage, together with an angular resolution 
of few milli-arcseconds. Simulations have shown that several 
hundreds of MBAs and TNOs wiU be observable with ALMA 
(see Busch 2009; Moufletet al. 2011). Interferometry at ther- 
mal wavelengths (mid-infrared to millimeter), in combination 
with other techniques {e.g., lightcurves, see Sect. 1.1), will thus 
allow the derivation of 3-D shape models, together with thermal 
properties, for many asteroids of small apparent diameter. 

2. The KOALA algorithm 

With advantages and drawbacks of each observing technique 
in mind, we have developped a multi-data inversion algorithm: 
Knitted Occultation, Adaptive-optics, and Lightcurve Analysis 
(KOALA), that makes simultaneous use of data from three dis- 
tinct observation types ^ to determine the physical properties 
of asteroids (Carry et al. 2010a; Kaasalainen 2011). KOALA 
takes advantage of the direct measure of the apparent size and 
shape of asteroids on the plane of the sky provided by the tim- 
ings of stellar occupations and disk-resolved images, and of the 
indirect constraints on spin and 3-D shape given by lightcurves. 
We quickly summarize below how the inversion works (see 
Kaasalainen 2011, for a comprehensive description of the al- 
gorithm). 



We categorize dense (i.e., lightcurves) and sparse ptiotometry togetlier to 
form a single data type. 
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KOALA is a step-iterative minimization algorithm, solv- 
ing for the spin parameters (spin-vector coordinates A, /3, 
and sidereal period P), 3-D shape (given by a set of N 
coefficients of spherical harmonics, including the overall 
size), phase function (defined as a three-parameter model, 
see Kaasalainen and Torppa 2001; Kaasalainen et al. 2001), 
and scattering law (generally taken as a combination of 
the Lommel-Seeliger and Lambert diffusion laws, following 
Kaasalainen and Torppa 2001, although other models such as 
Hapke can be used). 

From a set of estimated parameters (determined from 
lightcurve-only inversion or analysis of the disk-resolved im- 
ages, for instance), a trial solution is created, with associated 
synthetic data sets (i.e., simulated lightcurves, disk-resolved 
images, and occultation profiles). KOALA then follows a 
Levenberg-Marquardt minimization scheme to determine the 
set of parameters that best fit all the data sets simultaneously, 
by comparing at each step the synthetic data with real mea- 
surements. The iteration stops when the residuals between the 
simulated and measured data sets reach an acceptable level (i.e., 
the level of the intrinsic noise of the measurements). 

As a safeguard, the resolution (i.e., N) is set to the low- 
est possible value for which a fit to the data can be achieved. 
We also introduced several regularizations: (a) a non-convexity 
("smoothness") term, that avoids spurious features (unrealis- 
tic topography) at small scales, unconstrained by the data; and 
(b) an inertia tensor regularization for principal-axis rotators, 
that forces the asteroid spin-axis to remain aligned, within a 
few degrees, with the largest moment of inertia. The relative 
weights of the different data types and also of the regulariza- 
tions are determined using the maximum compatibility estimate 
of Kaasalainen (2011), instead of being subjective. 

3. KOALA and Rosetta-flyby shape model 

On 2010 July 10, the ESA Rosetta spacecraft made a close 
flyby of the main-belt asteroid (21) Lutetia. In support of 
the mission, we had combined optical lightcurves with disk- 
resolved images, from ground-based AO. We produced a full 
3-D shape model of Lutetia, using KOALA, months before the 
encounter (see Drummondet al. 2010; Carry et al. 2010b, and 
Fig. 1). This flyby provided a rare opportunity to test and cali- 
brate KOALA with close-up spacecraft imaging. 

The closest approach (CA) occurred at 3 170 km from the 
asteroid at a relative speed of 15km/s. The Narrow Angle 
Camera (NAC) of the OSIRIS instrument on-board Rosetta 
(Keller et al. 2007) returned a total of 202 images during the 
flyby, which spanned about 9 h. The NAC image scale ranged 
from 5 000 to 60 m/pix, reaching its minimum value at CA. 
These very high spatial-resolution images have been used to 
produce detailed 3-D shape models of Lutetia, using stereopho- 
toclinometry (Sierks et al. 2011), and stereophotogrammetry 
(Preusker et al. 201 1). We only consider here the model from 
Sierks et al. (2011), these two models being similar enough at 
the medium-to-large scale for our purpose. 

Lutetia's spin axis is tilted such that its pole is nearly 
in its orbital plane (obliquity of 96°, see Carry et al. 2010b; 







Diameters (km) 




Spin axis (°) 




d 


a b 


c 


A /3 ± cr 


KOALA^ 


105 


124+5 101+4 


93+13 


52 -6 +5 


OSIRIS ^ 


98 


121+1 101+1 


75+13 


52.2-7.8 + 0.4 



Table 1: Comparison of the dimensions and spin- vector coordinates of (21) 
Lutetia derived from OSIRIS images during the ESA Rosetta flyby (Sierks et al. 
2011) with those derived by KOALA (Drummond et al. 2010; Carry et al. 
2010b) prior to the flyby. Volume-equivalent diameter (d) and triaxial diam- 
eters (a > b > c) are reported in km, and spin-vector coordinates in degree 
(longitude A, latitude /? in the ECJ2000 reference frame, with cr the angular 
radius of the uncertainty circle). 



Sierks et al. 201 1). At the time of the Rosetta flyby, the south- 
ern hemisphere was in seasonal shadow, and observations at 
optical/near-infrared wavelengths were not possible south of - 
40° latitude. The detailed 3-D shape model derived from flyby 
images by stereophotoclinometry (Sierks et al. 2011) therefore 
does not cover a large fraction of the asteroid's southern hemi- 
sphere (see Fig. 2), the southernmost portion of the shape model 
being determined using the KOALA algorithm together with 
the ground-based data from Carry et al. (2010b). Almost all 
of the AO images that entered into the KOALA solution, even 
though from multiple epochs, were taken looking at either high 
southerly or high northerly sub-Earth latitudes (see Carry et al. 
2010b; Drummond et al. 2010). For both the ground-based 
KOALA and the flyby analyses, therefore, the shortest (c) di- 
mension (Table 1) is not well constrained and represents the 
largest contribution to the uncertainty in the volume estimates. 

We first present (Section 3.1) an overall comparison of 
our size and spin estimates (based on ground-based techniques) 
with those derived from the analysis of the very high spatial- 
resolution images acquired by OSIRIS (Sierks et al. 201 1). We 
then directly compare the predicted shape profiles, based on 
KOALA, with those images (Section 3.2), seeking to calibrate 
our method. 

3.1. Overall comparison 

In this section, we present a comparison of the spin-vector 
coordinates and dimensions derived from the flyby with the pre- 
flyby values from KOALA (Table 1). Because southerly lati- 
tudes were not visible during the flyby and the c-dimension in 
the OSIRIS model comes from a combination of OSIRIS imag- 
ing and KOALA model information, the OSIRIS c-dimension 
is not wholly independent. But the spin-vector coordinates and 
equatorial dimensions (a, b) from the OSIRIS analysis can be 
used as "ground-truth" to calibrate our KOALA method. 

First, the spin-vector coordinates agree to within two de- 
grees ^, well inside the 5° uncertainty quoted for KOALA Then, 
equatorial dimensions (a and b) are within 3 km of the flyby es- 
timates, again well within the uncertainties reported for each 
dimension using KOALA. The larger difference between the 
estimates of the short (c) axis results from the expansion of the 



'the coordinates obtained with KOALA are L2° and 1.8° from the spin- 
vector solutions derived by F. Preusker (A = 52.6°,/? = -7.1°, pers. communica- 
tion) and by Sierks et al. (2011) respectively, using the OSIRIS images. 
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Figure 2: Lutetia 3-D shape model from Sierks et al. (2011), displayed in a 
net layout. Regions imaged by OSIRIS are displayed in gray. The part of the 
shape model that is based on ground-based data only (lightcurves and images: 
LC+AO) is plotted in blue. 



KOALA c dimension from its original 80 km to 93 km, based 
on the analysis by Drummond et al. (2010). Indeed, the best- 
fit solution for all the AO images (including an additional ob- 
servation taken with a more equatorial geometry, presented by 
Drummond et al. (2010) and not used by Carry et al. (2010b) 
owing to calibration issues, see Carry et al. for details) pointed 
toward a larger c dimension (supported by independent con- 
siderations on the amplitude of lightcurves by Belskaya et al. 
2010). The volume-equivalent diameter is, therefore, larger 
(105 vs. 98 km) for the KOALA model than for the flyby- 
derived model by Sierks et al. (2011). Disk-resolved images 
of Lutetia taken with low sub-observer latitudes are required to 
confirm its c dimension, and to set tighter constraints on its vol- 
ume. 

The depth of large-scale concavities was slightly underesti- 
mated by KOALA (Fig. 3); estimating the depth of a concav- 
ity from profiles only is problematical because the concavity is 
hidden behind its rim. With stereophotoclinometry, is is possi- 
ble to sense depths. In addition, while KOALA is sensitive to 
large-scale concavities, the necessarily limited resolution when 
imaging from a distance of 200 million km means that it is less 
sensitive to small-to-medium-scale concavities. We evaluated 
the influence of craters on the volume of Lutetia. For that, 
we used the crater size distribution measured by Marchi et al. 




Figure 3: Similar plot as Fig. 2, showing the agreement between the pre-flyby 
KOALA (reddish, Carry et al. 201 Ob) and post-flyby (grayish, Sierks et al. 
201 1) shape models of (21) Lutetia. The radii of the KOALA model are larger 
in red regions, and those of the OSIRIS model in gray regions. The amalgam 
of gray and red colors illustrates the spectacular agreement between the re- 
sults from the Rosetta flyby and our KOALA model obtained before the flyby. 
The dominance of red hues, however, highlights the relativly lower sensitiv- 
ity to concavities of the KOALA model with respect to the OSIRIS 3-D shape. 
Shown in this comparison is the original KOALA model of Carry et al. (2010b), 
with a c-dimension of 80km (overall size 124x 101 X 80km), rather than the 
modified KOALA ("hybrid") estimate of Drummond et al. (2010), having a c- 
dimension of 93 km and shown in Table 1 . This is a more realistic comparison 
for the figure because much of the southern hemisphere of the OSIRIS model 
(Fig. 2) is based on the data that entered into the original KOALA model. 

(2011), and estimated their volume as that of a spherical cap 
(Vcrat), using the average depth-to-diameter ratio measured by 
Vincent et al. (2011): 

Vcr.t ^^d^^D^+d"^^ (1) 

where d is the depth of the crater and D its diameter We 
present in Fig. 4 the cumulative distribution of the volume en- 
compassed by craters, counted on four geomorphological units 
(see Thomas et al. 201 1, , for a detailed definition of the units), 
against their diameter. The vast majority of the volume is due 
to the handful of craters with diameters between 10 and 20 km. 
The total influence of these craters on the volume of Lutetia is 
0.6 + 0.1%, i.e., the volume of Lutetia would be 0.6% larger if 
these craters were not included in the 3-D shape model. Ex- 
trapolating this value, by considering the area covered by these 
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Figure 4: Cumulative distribution of the volume encompassed by craters as 
function of their diameter, for four geomorphological units of Lutetia: Achaia 
(Ac), Baetica (Bt), Narbonensis (Nb), and Noricum (Nr) (see the insert, adapted 
from Fig. 1 by Sierks et al. (201 1), and also Thomas et al. (201 1) for the defi- 
nition of these units). 

unit to the whole surface of Lutetia, make the total ^ influence 
of craters to be 2.4 + 0.6%. 

Proper modeling of the craters in the 3-D shape model (from 
high-resolution flyby images, Sierks et al. 2011) is therefore 
crucial, given the level of accuracy reached elsewhere on the 
surface. For Earth-based observations, however, this will re- 
main a minor source of uncertainty: depending on the methods, 
the volume accuracy ranges from few percent to few tens of 
percent (see Sect. 1). We discuss in the next section (3.2) the 
influence of these craters on the volume accuracy that can be 
reached using KOALA. 

3.2. Detailed analysis 

We push further the calibration of the KOALA method by 
comparing comprehensively the KOALA shape model predic- 
tions with the very high-spatial-resolution images provided by 
OSIRIS NAC. In the absence of a complete 3-D shape model, 
based on an entirely independent data set, it is difficult to fully 
calibrate the volume estimate provided by KOALA. For each 
OSIRIS image, we extracted the profile of the apparent disk of 
Lutetia, composed of its limb and its terminator We produced 
synthetic views of the KOALA shape model under the same 
geometry (i.e., as seen from Rosetta, see Fig. 5): phase angle, 
subsolar point (SSP) and sub-Rosetta point (SRP) coordinates, 
using the Miriade ^ VO ephemeris generator 

We estimate the fit of KOALA predictions to OSIRIS data 
by computing the difference between the projected areas on the 



Because the shape of the southern hemisphere is poorly constrained, and 
the number of craters there unknowfn, this extrapolated value is only a rough 
estimate. 

'http://vo.imcce.fr/webservices/miriade/ 

'"we used Rosetta flight kernel ORHR 00122.BSP 



2010-07-10713:48:18 
NAC 




Figure 5: Example of a profile comparison, as measured on the OSIRIS NAC 
detector plane (solid line) and simulated from the KOALA model (dotted line), 
taken at 13:48:18 UT (CA - lh57m). The light gray area represents the differ- 
ence in projected area on the plane of the sky between the prediction and the 
observation. We report the Rosetta-Lutetia distance, the coordinates of the sub- 
Rosetta point (SRP) and subsolar point (SSP), and a scale for angular (OSIRIS 
NAC pixels) and physical dimension. 

plane of the sky. The relative accuracy of the volume determi- 
nation (6V/V) can then be determined as follows: 

V 2 Jlo 

where and are the areas of the KOALA prediction and 
on OSIRIS frame, respectively. Negative and positive 6V re- 
spectively indicate an overestimate or underestimate of the vol- 
ume by the KOALA model. We discarded from the current 
analysis the images taken close to CA, which have substantial 
phase angle (above 10°), where local topography (e.g., crater 
rims) produce large projected shadows, increasing the uncer- 
tainty in the prediction of the terminator position on the surface 
of Lutetia. 

From this detailed analysis, we confirm the results from the 
overall comparison presented in Section 3.1: the volume was 
slightly overestimated with KOALA, relative to the results de- 
rived from the flyby. But, as visible in Fig. 6, the difference 
between the KOALA predictions and the OSIRIS images is al- 
most constant, at about -5%, for all images taken with phase 
angle smaller than ~10°. The largest deviations, still within the 
uncertainty reported using KOALA, are found for sub-Rosetta 
latitudes lower than about 40°, i.e., for geometries that had not 
been observed from Earth with AO. Owing to the restricted 
geometries of ground-based AO imaging observations to date 
(always close to "pole-on", see Drummond et al. (2010) and 
Carry et al. (2010b)), the differences between KOALA predic- 
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Figure 6: Difference in volume estimated (Eq. 2) from the difference in projected area between the profiles of Lutetia extracted from the OSIRIS images and the 
prediction from the KOALA model (Carry et al. 2010b). The symbols are color-coded as a function of the phase angle. Filled symbols correspond to images where 
the phase angle was smaller than 10°. 



tions and OSIRIS frames may be related to two distinct factors: 
one "inherent" to KOALA (which is the one we seek to eval- 
uate), and another related to the observing geometry. Unfortu- 
nately, there is no easy way to distinguish between these two 
effects. We thus estimate that the inherent uncertainty in the 
KOALA volume is about 5% (Fig. 6). This uncertainty can in- 
crease, however, due to unfavorable observing geometries, such 
as we have for the present case of Lutetia. 

In addition to calibrating the volume estimate, we assess the 
accuracy of the KOALA 3-D shape determination. From the di- 
rect comparison of both models (pre- and post-flyby, in Fig. 3), 
we can already qualitatively assert that KOALA allows accu- 
rate 3-D shape determination. However, many different shapes 
can result in similar volume, and also in similar overall triaxial- 
ellipsoid dimensions (which are accurate to a couple of km, see 
Section 3.1). We therefore apply the following criterion (6R, 
see Eq. 3) to estimate the local deviation of the KOALA model 
to the real shape of Lutetia, and therefore calibrate quantita- 
tively the KOALA 3-D shape determination: 

N ^ Oi ' 

where Oi and %, respectively, are the OSIRIS and KOALA 
profile radii of the ith point (out of describing the OSIRIS 
profile), measured from an arbitrary center 

We present in Fig. 7 the estimate of the local deviation of the 
KOALA predictions to the apparent shape of Lutetia measured 
on the OSIRIS NAC images, as a function of the time relative 
to CA. The deviations are limited to about 4 km at maximum, 
and are about 2 km on average (confirmed by the independent 
analysis of Preusker et al. 2011). The typical accuracy in the 
elevation of each vertex is therefore close to 2 km. This is con- 
sistent with the overall comparison presented in Section 3.1, 
but we show here that not only the overall sizes are accurate to 
about 2km, but that this accuracy is maintained at local scales. 



The apparent absence of craters in the KOALA model" 
(compared to 3-D shape models derived from flyby, e.g., 
Sierks et al. 2011; Preusker et al. 201 1) has therefore little in- 
fluence on the volume estimate and its accuracy. The largest 
craters seen on Lutetia have a diameter of about 20 km and 
a depth of about 2-3 km (see Fig. 4, and Marchi et al. 2011; 
Vincent et al. 2011), corresponding to the typical accuracy on 
elevation, and are therefore already included in the uncertainty 
envelop around the KOALA shape model. 

From this detailed comparison of OSIRIS frames with 
KOALA predictions, we assert that the KOALA model per- 
formed extraordinarily well, especially given two limiting 
factors related to the observing geometries of the ground-based 
observations. First, for the imaging observations, Lutetia 
had an angular diameter of about 0.1", close to the angular 
resolution for which we currently can extract useful shape 
information from ground-based telescopes with AO. Second, 
for the oppositions in 2007 and 2008, Lutetia was positioned 
at diametrically opposite apparent ecliptic coordinates from 
Earth, and because of the high obliquity, we observed Lutetia 
close to North-pole-on and then close to South-pole-on, and we 
were not able to achieve a good equatorial view, resulting in a 
poorly determined c-axis (Drummondet al. 2010; Carry et al. 
2010b). We expect to improve that in upcoming observations. 

We have determined the difference between the KOALA 
shape model and the real topography of Lutetia to be 2 km, 
on average. Considering Lutetia's volume-equivalent diameter 
(98 km, see Table 1), this translates into a relative precision on 
the radii of about 5%. We can therefore expect a conservative 
upper limit of 15% for the accuracy that can be achieved 
on volume estimates made by KOALA (including possible 
systematic effects), although it appears that the measured 
uncertainty (Fig. 6) is closer to 5%. These estimates of the 
accuracy achievable with KOALA are close to our previously 



"Although the model does not show round, crater- like, features, KOALA 
allows modeling of concavities, that can be either large impact craters. 
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Time relative to Closest Approach (h) 

Figure 7: Estimate of the local deviation (SR) between the profiles of Lutetia extracted from the OSIRIS images and the prediction from the KOALA model 
(Carry et al. 2010b). The symbols are color-coded as a function of the phase angle. Filled symbols correspond to images where the phase angle was smaller than 
10° (see text). 



estimated uncertainty: Although formal uncertainties for our 
shape-fitting algorithms are closer to 1 or 2 km, we had esti- 
mated that our size measurements were affected by systematic 
errors at the level of 1-3% (from simulations and observations 
of Saturn's moons, see Carry 2009; Drummond et al. 2009). 



4. KOALA and thermal radiometry 

As a supplementary investigation of the capabilities of 
KOALA for the study of a large sample of asteroids from 
ground-based observations, we analyze here the thermal prop- 
erties of (21) Lutetia, using only ground-based informa- 
tion: the KOALA 3-D shape model (Drummond et al. 2010; 
Carry et al. 2010b), and 104 individual thermal measurements 
(O'Rourke et al. 201 1, and reference therein), disregarding any 
information provided by the flyby of the asteroid by Rosetta. 
We then compare our results with those derived from the ther- 
mal infrared observations acquired by the MIRO instrument 
during Rosetta flyby (Gulkis et al. 2011). We also comment 
on the results obtained by O'Rourke et al. (201 1) with the same 
data set, but using the 3-D shape model derived from Rosetta 
imaging (Sierks et al. 201 1). 

Many observations at thermal wavelengths are available: 
Lutetia was observed multiple times by several infrared survey 
missions, like IRAS in 1983, or Akari in 2006/2007, but it was 
also targeted by different observing campaigns from ground: 
IRTF, ESO-T1MM12, and from space: Spitzer, Herschel (see 
O'Rourke et al. 2011, , for details on the observing circum- 
stances). 

First, we use the spin-vector coordinates provided by 
KOALA to determine the corresponding (lightcurve aver- 
aged) cross-sections for all observations that were used to de- 
rive the absolute magnitude of Lutetia: H = 7.25 + 0.01 mag 
(BoweUetal. 1989; Belskaya et al. 2010). All observations 
were taken with sub-observer latitude between -70° and -85°, 
close to pole-on geometry, when the apparent average diameter 



pv r (SI) References 

assumed 20 Gulkis et al. (201 1) 

0.185 + 0.005 5 O'Rourke et al. (2011) 

0.19 +0.01 <40 This work 



Table 2: Albedo (pv), and thermal ineilia (F, in SI units: Jm"^ s""^ K"') de- 
rived using MIRO on-board Rosetta (Gulkis et al. 201 1), ground-based data in 
combination with the shape model derived from Rosetta flyby (O'Rourke et al. 
201 1), and ground-based data in combination with the KOALA shape model 
(present study). 



of Lutetia was around 110-1 11 km (compared to its volume- 
equivalent diameter of 98 km, Sierks et al. 2011). The abso- 
lute magnitude of Lutetia was therefore slightly overestimated. 
We can still use it to determine the geometric visual albedo 
via the relation: \ogpv - 6.2559 - 21ogii - 0.4H 
(Pravec and Harris 2007, and references therein), providing d 
is the average apparent diameter of Lutetia at the time of the ob- 
servations. The combination of the KOALA shape model with 
the published absolute H magnitude leads to pv = 0. 19 + 0.01, 
in excellent agreement with the value derived during Rosetta 
flyby (0.19 + 0.01, see Sierks et al. 201 1). 

Second, the knowledge of the 3-D shape allows us to correct 
this bias on the actual H magnitude of Lutetia: Published ab- 
solute magnitudes are typically derived from a limited number 
of latitudes of the sub-observer point, and represent therefore 
an approximation only to the real absolute magnitude. Using 
the knowledge of the 3-D shape model with an absolute size 
scale and the geometric albedo (see above), one can determine 
the proper geometry-independent H-mag. Indeed, the absolute 
magnitude is intended to be a general, object-related property, 
rather than an observed quantity valid only for certain geome- 
tries (Bowelletal. 1989). The H-mag for Lutetia, based on 
the KOALA shape model and the geometric visual albedo of 
p-v = 0.19, is therefore H^ = 7.42 + 0.03. 

We use this refined H-mag to proceed with radiometric anal- 
ysis via a TPM code (Lagerros 1996, 1997, see Section 1.3). 



10 



This model considers a 1-D heat conduction into the surface, 
based on realistic surface conditions of illumination provided 
by the KOALA shape model (see O'Rourke et al. 2011, for 
details on such computations). As a complementary check of 
the techniques, we first determine the radiometric size and ge- 
ometric albedo of Lutetia, regardless of the size information 
provided by the KOALA shape model (following the method 
by Miiller et al. 2011). The calculation is based on the gen- 
eral thermal properties derived for large main-belt asteroids 
(Miiller and Lagerros 2002), a wavelength-dependent emissiv- 
ity model, a default thermal inertia F of 15 J nT^ s '^^ K and a 
default roughness implementation with 60% of the surface cov- 
ered by craters and an RMS of the surface slopes of 0.7 (see, 
e.g., Lagerros 1996, 1997, for definitions of these quantities). 
This leads to a radiometric volume-equivalent diameter (based 
on a dimensionless version of the KOALA 3-D shape model) 
of 99.8 +4.6 km and a geometric albedo pv of 0.198 + 0.017 
(weighted mean values and standard deviations from the anal- 
ysis of the 104 individual thermal measurements). These val- 
ues are in agreement with the pre-flyby estimates (diameter of 
98.3 ±5.9 km and geometric albedo pv of 0.208 + 0.025) by 
Mueller et al. (2006), using the lightcurve inversion model by 
Torppa et al. (2003). This demonstrates the value of combining 
thermal data with spin and 3-D shape models to derive accu- 
rately the diameters and albedos of asteroids. 

Finally, we constrain the thermal inertia of Lutetia by com- 
paring the TPM predictions (using the KOALA 3-D shape 
model, including its size estimate) for a range of thermal in- 
ertias with the 104 observed thermal measurements. This al- 
lows us to find the most probable thermal inertia to explain 
all data simultaneously, i.e., to match the before/after oppo- 
sition observations as well as the observations taken at the 
short-wavelength Wien part of the spectral energy distribu- 
tion (SED) and the long-wavelength Rayleigh-Jeans part of 
the SED. Thermal data also cover the entire rotation of Lute- 
tia and a significant range of sub-Earth latitudes (see the list- 
ing of observations in O'Rourke et al. 201 1). Thermal inertias 
in the range 5^0 J m"^ s "-^ K ' produce the lowest ;i'^-values 
(see Fig. 8), indicating the presence of fine-grained regolith 
(Brown and Matson 1987). These values are in agreement 
with the values measured by Rosetta (< 20 J m"^ s"°-^ K"', see 
Gulkis et al. 2011). The ;^'^-analysis also shows that one has 
to add substantial roughness on the surface to explain the ther- 
mal measurements. The KOALA 3-D shape model without any 
small-scale roughness cannot provide acceptable values. 

The prediction made using the combination of the KOALA 
and TPM models provides a good fit to all 104 measurements 
of Lutetia in the thermal infrared. There is a discrepancy, how- 
ever, for a given observing geometry that could be solved by the 
addition of a plateau/hill to the 3-D shape model, whose size re- 
mains within the quoted 3-D shape uncertainty as discussed by 
O'Rourke et al. (2011). 

From this detailed comparison of ground-based with flyby 
results, KOALA has been shown to be a powerful technique 
for the study of asteroids from the ground, not only for 
size/shape/spin, but now also for thermal properties. The val- 
idated (see Section 3) 3-D shape models, determined using 
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Figure 8: Thermal inertia determination using the KOALA shape model and 
ground-based observations. We use four values of surface roughness (default 
roughness means 60% of the surface is covered by craters and the RMS of 
the surface slopes is 0.7; see detail in O'Rourke et al. 2011) for thermal inertia 
ranging from very low to moderate (1-100 Jm"^ s""-^ K"'). 

KOALA, allow the study of the surface properties (e.g., albedo, 
thermal inertia) in great detail, and the results are consistent 
with those derived from the flyby. 

5. Concluding remarks: the future of KOALA 

The flyby of (21) Lutetia by ESA Rosetta provided a 
spectacular demonstration of the capabilities of the KOALA 
algorithm. Spin-vector coordinates were found to be accu- 
rate within two degrees, 3-D shape modeling to better than 
2 km (local topography), and dimensions to 2%. Volume 
estimates provided by KOALA are seen to be accurate to 
better than 10%. The thermal properties (albedo and thermal 
inertia) of Lutetia determined using a thermophysical model 
in conjunction with the KOALA shape model agree with the 
Rosetta-flyby-derived values, within the quoted uncertainties. 
These levels of accuracy on the spin and 3-D shape/size are 
typical for large main-belt asteroids, and not specific to Lutetia. 
Although it was extensively observed from the ground, being a 
spacecraft target, the number of lightcurves was similar to that 
of other large main-belt asteroids and the geometry of imaging 
observations was not particularly favorable (i.e., mostly close 
to "pole-on"). 

This ability of KOALA to determine the volume of main-belt 
asteroids of size ~ 100 km with an accuracy of about 10% opens 
the possibility for study of a larger set of small bodies using 
Earth-based observations. For instance, it can be expected 
to help efibrts to better understand the densities of asteroids 
belonging to different taxonomic classes (DeMeo et al. 2009). 
This will be assisted by mass estimates from about 150-200 
asteroids that will be determined from gravitational deflections 
observed by the upcoming ESA Gaia astrometry mission 
(Mouret et al. 2007), and the ever-growing number of known 
binary asteroids. During the first stages of the encounter 
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of Rosetta with its main target, comet 67P/Churyumov- 
Gerasimenko in 2014, we also plan to use a KOALA-style 
analysis of the first, low-resolution, resolved images to quickly 
produce a shape/size model. This will support the mission 
science until a full high-resolution model can be derived. 

The present implementation of KOALA allows the com- 
bined use of optical lightcurves (including sparse photometry), 
profiles from disk-resolved images, and chords from stellar 
occultations. We continue the development of KOALA toward 
the use of more data modes (e.g., interferometry, thermal 
radiometry) to increase the number of possible targets, and to 
set better constraints on targets observable only with certain 
techniques. Using different wavelength ranges, in particular 
in the thermal infrared, opens the possibility of deriving addi- 
tional physical properties, like albedo or thermal inertia. We 
also plan to incorporate additional cross-checks and constraints 
on the inversion, such as rigorous attention to differences in the 
spatial resolution of input images, to improve the confidence 
on the non-convex features and details. 
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